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\ Abstract. The formalism of the simple model of galactic chemical evo- 

. lution (GCE) and its main ingredients (stellar properties, initial mass 

(-H ' function, star formation rate and gas flows) are presented in this tuto- 
rial. It is stressed that GCE is not (yet) an astrophysical theory, but 

I ■ it merely provides a framework in which the large body of data con- 

O ' cerning the chemical composition of stars and gas in galaxies may be 

, intrepreted. A few examples illustrating those concepts are provided, 

^ ' through studies of the Solar neighborhood and of the Milky Way's halo. 

^ 1 Introduction 

> ; 

' Galactic Chemical Evolution (GCE) is the study of the evolution of the transfor- 
, mation of gas into stars and of the resulting evolution of the chemical composition 

■ of a galaxy. It has not yet the status of a full astrophysical theory, like e.g. the 
, theories of stellar evolution or of hierarchical structure formation in the universe. 
• The reason is our poor understandig of the driver of GCE, namely large scale star 

t — , formation in galaxies. On the contrary, the drivers of stellar evolution (energy pro- 

■ ducing nuclear reactions in plasmas heated by gravity) and of hierarchical structure 
^ , formation (cold dark matter assembled by gravity in an expanding universe) are 

! ■ fairly well understood. For that reason, those theories can make relatively robust 
rS \ predictions, concerning for instance, various evolutionary timescales, which is not 
' the case with GCE. 

Despite that, GCE provides a useful framework in which one may intrepret the 
large (and everexpanding) body of observational data concerning elemental and 
isotopic abundances in stars of various ages, as well as in the interstellar medium 
(ISM) and even the intergalactic medium (IGM). In particular, GCE allows one 
to: 

• Check/ constrain our understanding of stellar nucleosynthesis (expressed through 
the predicted yields of stars of various masses and metallicities) in a statisti- 
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cal way, i.e. by comparing GCE results to the mean trends and the dispersion 
of abundances and abundance ratios. 

• Establish a chronology of events in a given system, by finding when the 
metallicity reached a certain value, or when some stellar source (SNIa, AGBs 
etc.) became an important contributor to the abundance of a given nuclide. 

• Infer how a system was formed, by constraining the history of the SFR or 
of various gas movements, e.g. infall for the local disk or outflow for the 
Galactic halo (see Sec. 3). 

Different kinds of, progressively more sophisticated, models of GCE can be 
developed. In the simple model, ejecta of dying stars are instantaneously mixed 
in the ISM, which thus aquires a unique chemical composition at any time. Mod- 
els of inhomogeneous GCE relax this "instantaneous mixing approximation" in a 
semi-analytical manner and can account for dispersion of chemical abundances in 
stars and the ISM. Finally chemo- dynamical models treat also the dynamics of a 
galaxy's gas, stars and dark matter (either in a static or in a cosmological envi- 
ronment) providing a self-consistent framework for GCE (see e.g. Gibson et al. 
2003 for a review). One should keep in mind, however, that (i) both the number 
of free parameters and the number of potentially constraining observables increase 
with the degree of sophistication of the model, and (ii) all models, simple and 
sophisticated ones, suffer equally from our poor understanding of star formation 
(and stellar feedback in the case of chemo-dynamical models). 

In Sec. 2 we develop the formalism of the simple GCE model (Sec. 2.1) and 
we present the required ingredients: stellar properties (Sec. 2.2) and yields (2.3), 
stellar initial mass function (2.4), star formation rate (2.5) and gaseous flows (2.6); 
we also provide some useful analytical solutions of the system of GCE equations, 
obtained under the so called "Instantaneous recycling aprproximation" (Sec. 2.7). 
In Sec. 3 we illustrate the simple model by applying it to two well observed 
systems: the local disk (3.1) and the Milky Way's halo (3.2). 

2 Models of GCE: formalism and ingredients 

In the simple model of GCE, a galaxy consists initially of gas of primordial com- 
position (see contribution by A. Coc in this volume), namely by Xh '^0.75 for H 
and Xhe ~0.25 for '*He, as well as trace amounts of D, '^He and ^Li (abundances 
are given as mass fractions Xi for element or isotope i). The gas is progressively 
turned into stars with a Star Formation Rate (SFR) '^(t), with the star masses M 
having a distribution $(A/), called the Initial Mass Function (IMfJ^. Depending 



^In principle, the IMF may depend on time, either expHcitly or impHcitly (i.e. through a 
dependence on metallicity, which increases with time); in that case one should adopt a Star Cre- 
ation Function C{t,M) (making the solution of the GCE equations more difficult). In practice, 
however, observations indicate that the IMF does not vary with the environment, allowing to 
separate the variables t and M and adopt C{t,M) = ^{t)^{M). 
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on its lifetime tm, the star of mass M created at time t dies at time t + tm and 
returns a part of its mass to the interstellar medium (ISM), either through stellar 
winds (in the case of low mass and intermediate mass stars) or through supernova 
explosions (in the case of massive star^. The ejected material is enriched in el- 
ements synthesized by nuclear reactions in the stellar interiors, while some fragile 
isotopes (like D) are absent from its composition. Thus, the ISM is progressively 
enriched in elements heavier than H, while its D content is reduced. New stel- 
lar generations are formed from this ISM, their composition being progressively 
more enriched in heavy elements, i.e. with an everincreasing metallicity Z (where 
Z = 'EXi for all elements i heavier than He). 

In the framework of the simple model of GCE it is assumed that the stellar 
ejecta are immediately and efficiently mixed in the IS^M As a result, the ISM is 
caracterized at every moment by a unique chemical composition Xi{t), which is 
also the composition of the stars formed at time t. Since the surface composition 
of stars on the Main Sequence is not affected, in general, by nuclear reaction^ 
observations of stellar abundances reveal, in principle, the composition of the gas 
of the system at the time when those stars were formed. One may thus recover 
the chemical history of the system and confront observations to models of GCE. 



2.1 Formalism of the Simple Model of GCE 

The GCE scenario scetched in the previous paragraphs can be quantitavely de- 
scribed by a set of integrodifferential equations (see Tinsley 1980): 
The evolution of the total mass of the system m(t) is given by: 

If the system evolves without any input or loss of mass, the right hand member 
of Eq. (1) is equal to zero; this is the so called Closed Box Models the simplest 
model of GCE. The terms of the second member within brackets are optional and 
describe infall of extragalactic material at a rate f(t) or outflow of mass from the 
system at a rate o{t); both terms will be discussed in Sec. 2.6. 

The evolution of the mass of the gas mcit) of the system is given by: 

^ = -^ + E+[f-o] (2.2) 



■^Massive stars also eject part of their mass through a wind, either in the red giant stage (a 
rather neghgible fraction) or in the Wolf-Rayet stage (an important fraction of their mass, in the 
case of the most massive stars) . 

•^This is the so called Instantaneous Mixing Approximation, not to be confused with the 
Instantaneous Recycling Approximation (IRA), to be discussed in Sec. 2.7 

■^An exception to that rule is fragile D, already burned in the Pre-Main Sequence all over the 
star's mass; Li isotopes are also destroyed, and survive only in the thin convective envelopes of 
the hottest stars. 
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where ^{t) is the Star Formation Rate (SFR) and E{t) is the Rate of mass ejection 
by dying stars, given by: 

E{t) = (M - Cm) *(i - tm) $(M) dM (2.3) 

JMt 

where the star of mass M, created at the time t — tm, dies at time t (if tm < t) 
and leaves a Compact object (white dwarf, neutron star, black hole) of mass Cm, 
i.e. it ejects a mass M — Cm in the ISM. The integral in Eq. (3) is weighted by 
the Initial Mass Function of the stars $(M) and runs over all stars heavy enough 
to die at time t, i.e. the less massive of them has a mass Mt and a lifetime tm < t. 
The upper mass limit of the integral Mu is the upper mass limit of the IMF and 
is discussed in Sec. 2.4. 

Obviously, the mass of stars ms{t) of the system can be derived through: 

m = ms + Tog (2.4) 

As time goes on, a progressively larger part of the mass of stars is found in the 
form of compact objects of total mass mc{t), while the mass of stars still shining 
(live stars) is mL{t); one has then: 

ms = rriL + mc (2.5) 

The mass of compact objects is 



mc 



[ c{t) dt (2.6) 
Jo 



where c{t) is the mass creation rate of compact objects, given by: 



c{t) = / 

JMt 



Cm *(i - tm) $(M) dM (2.7) 



The mass of live stars mL{t) can then be calculated by using Eqs. 2.5 to 2.7. 

The evolution of the chemical composition of the system is decribed by equa- 
tions similar to Eqs. (2.2) and (2.3). The mass of element/isotope i in the gas is 
TOi = mcXi and its evolution is given by: 

^i^^^l^ = + E, + [fXij - (2.8) 

i.e. star formation at a rate ^ removes element i from the ISM at a rate ^Xi, 
while at the same time stars ro-injcct in the ISM that element at a rate Ei{t). If 
infall is assumed, the same element i is added to the system at a rate fXij, where 
Xij is the abundance of nuclide i in the infalling gas (usually, but not necessarily, 
assumed to be primordial). If outflow takes place, element i is removed from the 
system at a rate oXi^o where Xi^^ is the abundance in the outflowing gas; usually, 
o=Xj, i.e. the outflowing gas has the composition of the average ISM, but in 
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some cases it may be assumed that the hot supernova ejecta (rich in metals) leave 
preferentially the system, in which case Xj_o > for metals. 
The rate of ejection of element i by stars is given by: 



where Yi{M) is the stellar yield of element/isotope i, i.e. the mass ejected in 
the form of that element by the star of mass M. Note that Yi{M) may depend 
implicitly on time t, if it is metallicity dependent (see Sec. 2.3.2). 

The masses involved in the system of Eqs. 2.1 to 2.9 may be either physical 
masses, i.e. m, mc, ms etc. are expressed in and ^{t), E{t), c{t) etc. in 
M0 Gyr~^, or reduced masses {mass per unit final mass of the system), in which 
case m, mc, ms etc. have no dimensions and ^{t), E{t), c{t) etc. are in Gyr"-^. 
The latter possibility allows to perform calculations for a system of arbitrary mass 
and normalise the results to the known/assumed mass of that system; note that 
instead of mass, one may use volume or surface mass densities. 

The system of integrodifferential equations 2.1 to 2.9 can only be solved numer- 
ically, unless some specific assumptions are made, in particular the Instantaneous 
recycling aproximation (IRA, to be discussed in Sec. 2.7). Its solution requires 
three types of ingredients: 

• Stellar properties: stellar lifetimes tm, masses of stellar residues Cm and 

yields Yi{M)\ all those quantities can be derived from the theory of stellar 
evolution and nucleosynthesis and depend (to various degrees) on the initial 
stellar metallicity Z. 

• Collective stellar properties: the initial mass function $(M) and the Star 

Formation Rate none of them can be reliably derived from first principles 
at present, and one has to rely on empirical prescriptions. 

• Gas flows into and out of the system (infall, inflow, outfow, wind): in simple 
GCE models these factors are optional, i.e. their introduction depends on the 
nature of the considered galactic system (e.g. infall for the solar neighbor- 
hood or winds for small galaxies). In more physical (e.g. hydrodynamical) 
models, they should stem naturally from the physics of the system. 

We discuss those ingredients in the following sections. 
2.2 Stellar lifetimes and residues 

Stellar lifetime is a strongly decreasing fiuiction of stellar mass (see Fig. 1). Its 
precise value depends on the various assumptions (about e.g. mixing, mass loss, 
etc.) adopted in stellar evolution models and, most importantly on stellar metal- 
licity. Indeed, low metallicity stars have lower opacities and are more compact and 
hot than their high metallicity counterparts; as a result, their lifetimes are shorter 
(see Fig. 1 right). However, in stars with M>2 Mq, where H burns through the 
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Fig. 1. Left: Lifetimes of stars of solar metallicity Zq; points are from Geneva stellar 
models (Schaller et al. 1992) and the curve is a fit to those points (see Equ. 2.10). Right: 
Ratio of stellar lifetimes at Z©/20 and Zq for the same models. 



CNO cycles, this is compensated to some degree by the fact that the H-burning 
rate (proportional to the CNO content) is smaller, making the corresponding life- 
time longer; thus, for M>10 Mq, low metallicity stars live slightly longer than 
solar metallicity stars. Of course, these results depend strongly on other ingre- 
dients, like e.g. rotation (see Meynet, this volume). In principle, such variations 
in tm should be taken into account in GCE models; in practice, however, the 
errors introduced by ignoring them are smaller than the other uncertainties of the 
problem, related e.g. to stellar yields or to the IMfH. 

The lifetime of a star of mass M (in Mq) with metallicity Zq can be aproxi- 
mated by: 

T(Af) = 1.13 10^°A/"^ + 0.6 10^M"°-^5 + 1.2 10^ yr (2.10) 

This fitting formula is displayed as solid curve in Fig. 1 (left). A Zq star of 1 
Mq, hke the Sun, is bound to live for 11.4 Gyr, while a 0.8 Mq star for '^23 Gyr; 
the latter, however, if born with a metallicity Z<0.05 Zq, will live for "only" 13.8 
Gyr, i.e. its lifetime is comparable to the age of the Universe. Stars of mass 0.8 
Mq are thus the lowest mass stars that have ever died since the dawn of time (and 
the heaviest stars surviving in the oldest globular clusters). 

The masses of stellar residues are derived from stellar evolution calculations, 
confronted to observational constraints. In the regime of Low and Intermediate 
Mass Stars (LIM^, i.e. for M/M0<8-9, the evolutionary outcome is a white 



^Metallicity dependent lifetimes have to be taken into account in models of the spectro- 
photometric evolution of galaxies, where they have a bigger impact. 

^LIMS are defined as those stars evolving to white dwarfs. However, there is no universal 
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Fig. 2. Left: Masses of stellar residues as a function of initial stellar mass, for stars 
of metallicity Zq; for massive stars (M>30 M0) the two curves correspond to different 
assumptions about mass loss, adopted in Limongi and ChiefR {solid curve, this volume) 
and Woosley and Ifeger (2007, dotted curve), respectively. Right: Mass fraction of the 
ejecta as a function of initial stellar mass; the two curves for M>30 M© result from the 
references in the left figure. 



dwarf (WD), the mass of which (in M0) is given by (Weidemann 2000): 

Cm{WD) = 0.08 M + 0.47 {M < 9) (2.11) 

Thus, the heaviest WD has a mass of ^1.2 M0, slightly below the Chandrasekhar 
mass limit of ^^1.44 Mq. The less massive WDs formed from single mass stars 
have a mass ^^0.53 Mq. Observations show indeed that WD masses cluster around 
that value. Note that WD with much lower masses have been found (even below 
0.2 Mq, e.g. Kilic et al. 2006), but they result from binary system evolution. 

Stars more massive than 8-9 Mq explode as supernovae (SN), cither after 
electron captures in their 0-Ne-Mg core (M<11 Mq) or after Fe core collapse 
(M>11 Mq). The nature and mass of the residue depends on the initial mass 
of the star and on the mass left to it prior to the explosion. It is often claimed 
that solar metallicity stars of M<25 Mq leave behind a Neutron Star (NS), while 
heavier stars leave a black hole (BH). Neutron star masses are well constrained by 
the observed masses of pulsars in binary systems: Mjvs=1.35±0.04 Mq (Thorsett 



definition for the mass limits characterizing Low and Intermediate Mass stars. The upper limit 
is usualy taken around 8-9 Mq , although values as low as 6 Mq have been suggested (in models 
with very large convective cores). The limit between Low and Intermediate masses is the one 
separating stars powered on the Main Sequence by the p-p chains from those powered by the 
CNO cycle and is ~1.2-1.7 Mq, depending on metallicity. 
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and Chakrabarty 1999), which is adopted here 

Cm{NS) = 1.35 (9 < M < 25) (2.12) 

i.e. Cm is independent of the initial mass M in that case. However, black hole 
masses are not known observationally as a function of the progenitor mass, while 
theoretical models are quite uncertain in that respect. Thus, for solar metallicity 
stars with mass loss, the Frascati models (Limongi, this volume) produce black 
hole masses proportional to the initial mass, which can be approximated (in Mq) 
by: 

Cm{BH) = 0.24 Af -4 (M>25) (2.13) 

while Woosley and Hegcr (2007) find that, as a result of high mass losses for M>25 
M0, the average black hole mass (in M©) is 

Cm{BH) =3 (Af > 25) (2.14) 

(see Fig. 2). It is commonly accepted that black hole masses are much larger at 
low metallicities, where the effects of mass loss are negligible. Note, however, that 
the magnitude of the effect could be moderated in models with rotational mixing, 
which induces mass loss even at very low metallicities (see Meynet, this volume). 

It should be stressed that the mass limit between the progenitors of neutron 
stars and black holes is not understood at present, and adopting simple theoretical 
"recipes" can be strongly misleading. For instance, Muno (2006) points out that 
a magnetar (strongly magnetised neutron star) is found in a star cluster less than 
4.5 Myr old (as deduced from the absence of red supergiants in it), implying a 
progenitor mass larger than 50 Mq. Of course, it is possible that strong mass loss 
in a binary system reduced the mass of such a heavy progenitor to a very low level 
prior to the explosion. 



2.3 Stellar yields 
2.3.1 Definitions 

The quantities required in Eq. 2.9 are the stellar yields Yi{M), representing the 
mass ejected in the form of element i by a star of mass M . Those quantities are 
obviously Yi{M) >0 (Fj=0 in the case of an isotope totally destroyed in stellar 
interiors, e.g. deuterium). However, their value is of little help in judging whether 
star Af is an important producer of isotope i (e.g. by knowing that a 20 Mq star 
produces 10""^ Mq of Mg or 1 Mq of O, one cannot judge whether such a star 
contributes significantly - if at all - to the galactic enrichmemnt in those elements). 

More insight in that respect is obtained through the net yields yi{M), which 
represent the — it newly created mass of nuclide i from a star, i.e. 

y,{M) = Y,{M)- ModM) (2-15) 

where Mo^i{M) is the mass of nuclide i originally present in the part of the star 
that is finally ejected: 



Mo,i{M) = Xo,i{M - Cm) 



(2.16) 
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Table 1. Yield definitions 



Nuclide i 


Yields Y,{M) 


Net yields yi{M) 


Production factors fi{M) 


Created 




> 


> 1 


Re-ejected 


= Mo,i 


= 


= 1 


Destroyed 


< Mo,, 


< 


< 1 



Mo^i is defined in Eq. 2.16. 



and Xo^i is the mass fraction of nuclide i in the gas from which the star is formed. 
Obviously, yi{M) may be positive, zero or negative, depending on whether star M 
creates, simply re-ejects or destroys nuclide i. Net yields are not usedin numerical 
models of GCE, but they are usedin analytical models, adopting the Instantaneous 
recycling approximation (see Sec. 2.7). 

Finally, the production factors fi{M) are defined as: 

They are useful in the sense that they immediately indicate whether star M is 
an important producer of nuclide i. For instance, massive stars are the exclusive 
producers of oxygen, for which / '^lO on average (see Fig. 4). If such stars produce 
another nuclide L with, say, / ^3 only, they are certainly important contributors, 
but they cannot account for the solar L/O ratio; another source is then requireqj 
for nuclide L. Note that the use of production factors as defined in Eq. 2.17 and 
described in the previous paragraph is interesting only when comparison is made 
for a star of a given initially metallicity (see e.g. Fig. 4, where denominators 
in Eq. 2.17 are calculated for Zq always, even if yields are for different initial 
metallicities) . 

The properties of the various quantities defined in this section are summarized 
in Table 1. 



2.3.2 Metallicity effects: primary vs. secondary and odd vs. even nuclides 

Nuclides with yields independent of the initial metallicity of the star are called 
primary. They are produced exclusively from the initial H and He of the star (the 
abundances of which vary little throughout galactic history). This is, in principle, 
the case for most isotopes up to the Fe peak, as well as the r-isotopes. In stellar 
models with no mass loss, such as those constructed until recently (including the 
widely used Woosley and Weaver 1995 models) initial metallicity plays a negligible 
role in the overall structure of the star, and thus the yields of primary elements 
(like e.g. oxygen or calcium) correspond to their definition, i.e. they are metallicity 
independent. However, in models with mass loss the situation is different. Nuclides 



^The example is taken on the case of iron, for which another source is required, beyond 
massive stars; that source is SNIa (see Sec. 2.3.5) 
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produced early in the lifetime of the star (i.e. during H and He burning), like C and 
O, are directly affected by mass loss: for instance, large amounts of material are 
ejected in the form of He and C in the Wolf-Rayet phase of high metallicity stars, 
leaving a small core available to turn into O; on the contrary, lower metallicity stars 
produce more O and less He and C (see Meynet, this volume). Thus, the yields 
of those nuclides, although defined as primaries, show a considerable dependence 
on metallicity. Nuclides produced in later evolutionary phases, much closer to the 
stellar center (like Si and Ca) are less affected by mass loss. Note that rotational 
mixing, inducing mass loss even at low metallicities, makes the mass loss effect 
less metallicity dependent. 

Nuclides with yields dependent on the abundance of some non-primordial nu- 
cleus in the original stellar composition are called secondary. This is (again, in 
principle !) the case of : s-nuclides, produced by neutron captures on initial Fe- 
peak nuclei ; ^''N and ^^C in moderately massive stars (made from initial C and 
O in the CNO cycle of H-burning); ^^O, made by ^^N (a, ^^^F {13+^^0 in early 
He-burning inside massive stars. Note that in the latter case '^^O is made not from 
the initial ^''N of the star (itself a secondary nuclide), because in that case ^^O 
would be a tertiary nuclide; it is made from the ^^N already synthesised in the star 
through the CNO cycle, i.e. it is essentially made from initial C and O. In fact, 
we know of no tertiary nuclide, made from a secondary one initially present in the 
star. For instance, p-nuclides are made from the photo-desintegration of s-nuclei 
(see Arnould, this volume) which are produced inside the star from the initial Fe 
(a primary nuclide); thus, p-nuclides are also secondary. 

Most of the metals in a massive star (namely C and O) turn into ^^N in H- 
burning and most of that "'^^N turns into ^^Ne during He-burning, through the 
reaction chain leading to (see previous paragraph) plus ^*0(a,7)^^Ne. The 
beta-decay involved in this chain modifies the neutron-to-proton ratio of the star, 
since the initial mixture of -"^^C and "^^O has the same number of neutrons N and 
protons Z (its neutron excess rj = is 77=0) while ^^Ne (Z=10 and N=12) has 
77=0.09. This, metallicity dependent, surplus of neutrons affects the products of 
subsequent burning stages: odd nuclides like e.g. ^''Na, ^'''Mg, ^^Al etc. are more 
readily produced in large rj (large initial metallicity) environments with respect to 
even nuclides; this the so-called odd-even effect (Arnett 1996). For such nuclides, 
use of metallicity- dependent yields is mandatory for GCE studies. 

2.3.3 Yields of massive stars 

Massive stars produce the quasi-totality of nuclides between carbon and the Fe- 
peak, as well as most of the heavier than iron nuclei: the light s- (up to Y), the r- 
(neutron-rich) and p- (neutron poor, w.r.t. to the nuclear stability valley) nuclides. 

The typical composition profile of a massive star after its supernova explosion 
is shown in Fig. 3. It displays the typical onion-skin structure, with heavier 
elements encountered in its various layers from the surface to the centc;r. The 
thickness of the layers depends on assumptions about convection and other mixing 
processes. The abundances of the nuclides in each layer depend essentially on 
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Fig. 3. Interior composition of a 25 Mq star after its explosion; only major isotopes are 
displayed (from Woosley and Weaver 1995). 



adopted nuclear reaction rates (but also on the mixing processes). The abundances 
in the outer layers depend also on the adopted mass loss rates (combined to the 
mixing processes). The abundances in the inner layers depend also on assumptions 
about the explosion of the star (mass cut, energy, fallback etc.). In particular, 
intermediate mass elements produced in massive stars may be divided in three 
main classes: 

- In the 1st class belong N, C, O, Ne and Mg, which are mainly produced in 
hydrostatic burning phases and are found mainly in layers that are not heavily 
processed by explosive nucleosynthesis; the yields of those elements depend on 
the pre-supernova model (convection criterion, mixing processes, mass loss and 
nuclear reaction rates). In particular, rotational mixing may produce primary N, 
even in very low metallicity stars (see Meynet, this volume). 

- In the 2nd class belong Al, Si, S, Ar and Ca. They are produced by hydrostatic 
burning, but their abundances are substantially affected by the passage of the 
shock wave. Their yields depend on both the pre-supernova model and the shock 
wave energy. 

- In the 3cl class belong the Fe-peak nuclei, as well as some lighter elements like 
Ti; their yields depend crucially upon the explosion mechanism and the position 
of the "mass-cut" . 
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Fig. 4. Production factors fi{M) = ^ ^'j^j^^g^^) for elements up to Ge in the ejecta 
of stars with no mass loss, of initial metallicity Z/Zq=1, 0.1, 0.01 and 0.0001 (upper 
left panel) in the mass range 12-35 Mq (data from ChiefS and Limongi 2004). Primary 
elements (C, O, Ne, Mg, Si, S, Ar, Ca, Ti, Cr, Fe) have metallicity independent /. N is 
produced only as secondary in massive stars with no rotational mixing. The odd-even 
effect is clearly seen for Na, Al, P, CI, K, Sc, V, Mn,while Li, Be and B are depleted 
to various degrees. A horizontal dotted line at f=l in all panels indicates ejection of 
material with solar abundance. 



In view of the important uncertainties still affecting the ingredients of stellar 
physics, it is obvious that stellar yields have large uncertainties. How could the 
validity of the theoretical stellar yields be checked? Ideally, individual yields should 
be compared to abundances measured in supernova remnants of stars with known 
initial mass and metallicity! However, such opportunities are extremely rare. In 
the case of SN1987A, theoretical predictions for a 20 M0 progenitor are in rather 
good agreement with observations of C, O, Si, CI and Ar (Thielemann et al. 
1996). SN1987A allowed also to "cahbrate" the Fe yield (~0.07 Mg) from the 
optical light curve (powered at late times from the decay of ^^Co, the progeny of 
^^Ni), extrapolated to the moment of the explosion (e.g. Arnett et al. 1989). This 
may be the best way to evaluate the Fe yields of other SNII at present, until a 
convincing way of determining the "mass-cut" from first principles emerges. 
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Our confidence to massive star yields stems rather from the fact that, when 
averaged over a reasonable IMF and used in a GCE model with a SFR appropriate 
for the solar neighborhood, they reproduce the solar composition between C and 
Fe-peak within a factor of 2-3 (see Sec. 3.1.2). This is a remarkable achievement, 
taking into account that solar abundances in that mass range vary by factor of 
~10^ between the most abundant of them (O) and the less abundant one (Sc). A 
similar situation holds for the light s-nuclides and the p-nuclides (both products 
of massive stars), whereas the situation with r-nuclide yields is problematic at 
present (see Arnould, this volume). 



2.3.4 Yields of intermediate mass stars 

Intermediate mass stars synthesize substantial amounts of several nuclides, mainly 
in the Asymptotic Giant Branch (AGB) phase, when H and He burn intermittently 
in two shells surrounding the inert CO core. Nucleosynthesis occurs in those 
shells, as well as in the bottom of the convective AGB envelope, if it penetrates 
in regions of high enough temperature {Hot Bottom Burning or HBB). Such stars 
are the main producers of heavy s-nuclei at a galactic level and they synthesize 
large amounts of ^He, ^^N, ^^C, ^''O, ^^F. etc. (see Siess, this volume). In fact, 
they were thought to be the main producers of those light nuclides. However, 
massive stars can also produce large amounts of those nuclides (especially when 
rotational mixing is taken into account) and eject them through stellar winds (see 
Meynet, this volume); in view of the uncertainties affecting the yields of both 
AGBs and massive stars, it is not clear yet whether the contribution of AGBs to 
the solar abundances of those nuclides is predominant or not. For instance, it has 
been suggested that massive stars are the main source of carbon in the Milky Way 
(Prantzos et al. 1994). 

Yield uncertainties are, in general, more important in the case of AGBs than 

in the case of massive stars. The role of mixing is indeed critical, in several 
respects: for providing the appropriate neutron source in the case of the s-process 
(i.e. through ^^C(Q!,n) in the intershell region, see e.g. Goriely and Siess 2006); for 
bringing material from the He-shell up to the envelope (the so-called third dredge- 
up); and for bringing the bottom of the convective envelope to sufficiently high 
temperatures as to burn hydrogen through the CNO cycle (HBB, synthesizing 
primary "'^''N and ^'^C). Mass loss is also critical, since it determines the duration, 
and the degree of nuclear processing, of the AGB phase. Note that, due to the 
complexity of the AGB phase, yields are often calculated by "post-processing", 
i.e. nucleosynthesis is calculated in an envelope decoupled from the stellar core. 

Some results of nucleosynthesis calculations for AGBs are displayed in Fig. 5 
(integrated over an IMF) as a function of initial stellar metallicity. It is clearly seen 
that ^^N and ^"'C are produced as primaries (because of Hot Bottom Burning), 
whereas ^^O is depleted (net yields are negative). Yields of s-nuclides in AGBs as 
a function of metallicity are given in Goriely and Mowlavi (2000). 
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Fig. 5. Net yields of typical products of intermediate mass stars, integrated over an IMF 
and displayed as a function of initial stellar metallicity; data from various authors are 
compiled by Marigo (2001). 



2.3.5 Yields of thermonuclear supernovae (SNIa) 

Thermonuclear supernovae are exploding white dwarfs in binary systems (see Isern, 
this vohimc). Although neither the progenitor stars nor the exact mode of the 
propagation of the nuclear flame are well understood at present, it is generally 
believed that a large fraction of the white dwarf mass burns in conditions of Nuclear 
Statistical Equilibrium (NSE). In those conditions, large amounts of Fe-peak nuclei 
and, in particular, ^^Ni, are produced: about 0.6 Mq of ^^Ni, with extreme values 
ranging from up to 1 Mq down to 0.1 Mq. 

In most models, only the internal parts of the white dwarf burn all the way to 
NSE; intermediate mass nuclei (from Si to Ca) are also produced in the explosion, 
in agreement with observations of early spectra of SNIa. The most widely used 
model is the deflagration model W7, which reproduces satisfactorily the spectro- 
scopic requirements. The results of a recent version of that model (Iwamoto et al. 
1999) appear in Fig. 6. One may note a substantial production of Si, S, Ar and Ca 
(at about 30% of the corresponding Fe production) and a slight oveproduction of 
^"^Fe and ^®Ni (the major default of the W7 model on grounds of nucleosynthesis). 

SNIa are major contributors to the galactic content of ^^Fe, the stable product 
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Fig. 6. Production factors / for isotopes from C to Cu in the "canonical" model W7 of 
SNIa (from Iwamoto et al. 1999). They are normalized to /Fe56 = l- 



of ^^Ni decay, and in Fe-peak nuclides in general. In the case of the Milky Way, 
this can be seen as follows: The observed frequency of SNIa in external galaxies 
of the same morphological type (i.e. Sbc/d) is about 5 times smaller than the 
corresponding frequency of core collapse SN (SNII+SNIb,c), see e.g. Mannucci et 
al. (2005). But core collapse SN produce, on average '^O.l Mq of ^^Ni, that is 6 
times less than SNIa. Thus, SNIa contribute at least as much as massive stars to 
the production of Fe and Fe-peak nuclides in the Milky Way. In the case of the 
solar neighbor hod, this is coroborated by another observational argument, namely 
the evolution of the 0/Fe ratio, that will be presented in Sec. 3.1.2. 

SNIa may have long-lived progenitors, with lifetimes of up to several Gyr; this 
introduces a substantial delay in their rate of element ejection in the ISM. The evo- 
lution of SNIa rate depends on the assumptions made about the progenitor system 
and it obviously cannot be simply proportional to the SFR. One of the simplest 
parametrizations, introduced recently, consists in assuming that the SNIa rate is 
proportional to both the SFR and the cumulative mass of the stellar population 

RsNiait) = a^{t) +(3ms{t) (2.18) 

where ^(t) and ms{t) are defined in Eq. (2.2) and (2.4), respectively. The values 
of a and /? are adjusted, as to account for the observations of SN rates in external 
galaxies of different morphological types. It should be noted, however, that other 
formulations provide an equally satisfactory fit to the data. This is the case of the 
widely used parametrisation of Greggio and Renzini (1983), taking into account 
the time delay introduced by the evolution of the secondary star and the mass 
distribution of the two components of the binary system. Despite their uncer- 
tainties, due to our poor understanding of the progenitor systems of SNIa, such 
formulations offer more physical insight into the problem. 
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2.4 Initial Mass Function (IMF) 

The IMF cannot be calculated at present from first principles (despite a large body 
of theoretical work, see e.g. Bonnell et al. 2006 for a recent review) and it has to 
be derived from observations. However, this derivation is not quite straightforward 
and important uncertainties still remain, especially in the region of massive stars. 

In a first step, from the observed Luminosity function f{L) = dN/dL and the 
theoretical mass-luminosity relation for Main sequence stars <?i>(i) = dM/ dL , the 
Present Day Mass Function F{M) is derived 

, , dN dL , , 

F(M) = — — 2.19 
^ ^ dL dM ^ ' 

Assuming that (l){L) is well known, one has still to account for various biases in the 
derivation of F{M), like e.g. the role of binary and multiple stellar systems. In the 
case of stellar clusters, the impact of cluster dynamics (leading to the preferential 
evaporation of the lowest mass members) has to be taken into account. In the case 
of the field F{M) in the solar neighborhood, one has to account for the fact that 
the dispersion of stellar velocities perpendicularly to the galactic plane increases 
with time, i.e. the corresponding scale- height of the stellar population h(L) is a 
decreasing function of L. 

Once the bias-corrected PDMF is established, the IMF may be derived, in 
principle, if the Star formation rate (SFR) is known. In fact, the PDMF is the 
integral over time of the star creation rate C{M, t) = $(M)\l/(t), from which dead 
stars are removed. Assuming that the IMF deos not depend on time, one has: 

m) = (2-20) 

where T is the age of the system (T ~12-13 Gyr in the case of the solar neighbor- 
hood and the Milky Way disk). That expression can be simplified in two cases, 
where knowledge of ^'(t) is not really required: 



In the case of short-lived stars (tm «T, i.e. M > 2 Mq) 



HM) = (2.21) 

where is the current SFR, assumed to be constant during the last Gyr. 

• In the case of "eternal" stars (tm >> T, i.e. M < 0.8 Mq for low metallicity 
stars) 

$(M) = f f ^ (2.22) 

where < * >= ^{t)dt/T is the past average SFR during the system's 
history. 
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Fig. 7. Slope 1 + X of the IMF (assumed to be described by a multi-power-law form), 
according to observations in various asrophysical environments; the dashed horizontal 
lines indicate average values in three selected mass ranges, with 1 -|- X=2.35 being the 
classical Salpeter value (from Kroupa 2002). 



Note that the shape of the IMF $(M) in Eq. (2.21) or (2.22) does not depend on 
the values of 5*0 or < vj/ >. For the intermediate mass regime (0.8< M/Mq <2) 
the situation is more complicated and knowledge of \l/(i) is, in principle, required. 
Continuity arguments may be used to bridge the gap between the low and high 
mass regimes. The evolution of the PDMF in the solar neighborhood is illustrated 
in Fig. 8 (left panel), where it is also compared to field star observations (right 
panel) . 

Based on observations of stars in the solar neighborhood and accounting for 
various biases (but not for stellar multiplicity), Salpeter (1955) derived the local 
IMF in the mass range 0.3-10 Mq, a power-law IMF: 

$(M) = 4^ = AAf-(i+^) (2.23) 

with slope Jil=1.35. That slope is indeed found in a large variety of conditions 
and the Salpeter IMF is often used in the whole stellar mass range, from 0.1 to 
100 Mq, especially in studies of the photometric evolution of galaxies. However, 
it is clear now that there are less stars in the low mass range (below 0.5 M©) than 
predicted by the Salpeter slope of X=1.35. For instance, Kroupa (2002 and Fig. 
7) adopts a multi-slope power-law IMF with X=0.35 in the range 0.08 to 0.5 M©, 
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Fig. 8. Left: Evolution of the Present Day Mass Function, displayed here at four different 
moments (0, 1, 4 and 13 Gyr) for a smooth star formation history, correspondig to the 
one of the solar neighborhood. The initial PDMF (0 Gyr) is the same as the adopted 
IMF, while the final PDMF (13 Gyr) is the one observed today. The latter curve is 
compared on the right panel with the PDMF derived for the solar neighborhood by Scalo 
(1986, error bars). The adopted IMF is log-normal below 1 M© and a power-law with 
slope above 1 Mq. 



whereas Chabrier (2003) prefers a log- normal IMF below 1 Mq. 

For GCE purposes, low mass stars are "eternal" and just block matter from 
recycling in the ISM. Most important, in that respect, is the shape of the IMF 
above 1 M0. Unfortunately, the situation is not clear yet. Observations of the 
IMF in various environments, and in particular young clusters (where dynamical 
effects are negligible) suggest that a Salpeter slope X=1.35 describes the data 
well (see also Fig. 7). However, determination of the IMF in young clusters 
suffers from considerable biases introduced by stellar multiplicity and pre-main 
sequence evolution. Moreover, in the case of the field star IMF in the the solar 
neighborhood, Scalo (1986) finds X=1.7, i.e. a much steeper IMF than Salpeter. 

According to Weidner and Kroupa (2006), the slope of the stellar IMF is indeed 
the one observed in clusters, but GCE studies involve the galaxian IMF, i.e. the 
sum of all cluster IMFs, which is steeper than the stellar IMF. The reason is that, 
although every single cluster has the same stellar IMF (X=1.35), the maximum 
stellar mass Mmax,c in a a cluster increases with the total mass of that cluster: 
indeed, observations reveal that small clusters may have Mmax.c as low as a 
few Mq, whereas large clusters have Mmax,c up to 150 Mq. If this were just 
a statistical effect, the slope of the resulting galaxian IMF (i.e. the sum over 
all cluster IMFs) would also be X=1.35, but if there is a physical reason for the 
observed Mmax,c vs Mciuster relation, then the resulting galaxian IMF would 
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Fig. 9. Top: Three initial mass functions: solid curve: Salpeter (power-law in the whole 
mass range), dotted curve: Kroupa (multi-slope power law for M<1 M0) -I- Scalo {X—1.7 
for M>1 M0), dashed curve: Chabrier (log-normal for M<1 Mq) -|- Salpeter (X=1.35 
for M>1 M0). Bottom: Ratio of the three IMFs to the one of Salpeter. 



necessarily be steeper. In that way, one may conciliate the results of Scalo (1986, 
a study still unique in its kind) with observations of cluster IMFs, which favour 
the Salpeter valutH. For GCE studies, the Scalo value of X=1.7 appears more 
appropriate, and it is indeed used in most detailed models of the local GCE. 
The IMF is normalised to 

/•Mu 

$(M) = / $(Af) M dM ^ 1 (2.24) 

J Ml 

where Mu is the upper mass limit and Ml the lower mass limit. Typical values are 
Mjj ~100 Mq and Ml ~0.1 M0, and the results depend little on the exact values 
(if they remain in the vicinity of the typical ones). A comparison between three 
normalized IMFs, namely the "reference" one of Salpeter, one proposed by Kroupa 
(with the Scalo slope at high masses) and one by Chabrier (with the Salpeter slope 
at high masses) is made in Fig. 9. 



*Some implications of that idea, suggesting that the slope of the high mass part of the IMF 
is variable and depends on the star forming activity of the galaxy, are explored in Koppen et al. 
(2006). 
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A useful quantity is the Return Mass Fraction R 

r-Mu 

R = (M - Cm) HM) dM (2.25) 

J Mt 

i.e. the fraction of the mass of a stellar generation that returns to the ISM. For the 
three IMFs displayed in Fig. 9 one has R= 0.28 (Salpeter), 0.30 (Kroupa+Scalo) 
and 0.34 (Chabrier+Salpeter) , respectively, i.e. about 30% of the mass gone into 
stars returns to the ISM. 

2.5 Star Formation Rate 

Star formation is the main driver of galactic evolution and the most uncertain pa- 
rameter in GCE studies. Despite decades of intense observational and theoretical 
investigation (see e.g. Elmcgreen 2002 and references therein) our understanding 
of the subject remains frustratingly poor. 

Observations of various SFR tracers in galaxies provide only relative values, 
under the assumption that the IMF is the same everywhere. Most such tracers 
concern formation of stars more massive than ~2 M.q ; very little information exists 
for the SFR of low mass stars, even in the Milky Way. 

Some commonly used tracers of SFR (i.e. of young stellar populations) are: 

• Luminous star counting in some region of the Hcrtzsprung- Russell diagram; 
this method can be applied only to nearby galaxies like the Magellanic clouds. 

• Flux measurements in recombinations lines H^, etc., emitted by gas 
which is ionized by nearby OB stars (emitting photons bluewards of 912 A). 

• UV flux of OB stars, which are relatively young (the UV emissivity of Hori- 
zontal branch stars and of planetary nebulae nuclei is negligible, in general). 

• Far infrared (FIR) emission of dust, heated by absorption of UV photons 
emitted by hot nearby stars. 

• Surface density of supernova remnants, like pulsars (in star forming galaxies, 
i.e spirals and irregulars). 

Those tracers reveal that star formation occurs in different ways, depending on 
the type of the galaxy. In spiral disks, star formation occurs mostly inside spiral 
arms, in a sporadic way. In dwarf, gas rich, galaxies, it occurs in a small number 
of bursts, separated by long intervals of inactivity. Luminous Infrared galaxies 
(LIRGS) and starbursts (as well as, most probably, ellipticals in their youth) are 
caracteriscd by an intense burst of star formation, induced by the interaction (or 
merging) with another galaxy. 

Today there is no theory to predict large scale star formation in a galaxy, 
given the various physical ingredients that may affect the SFR (e.g. density and 
mass of gas and stars, temperature and composition of gas, magnetic fields and 
frequency of collision between giant molecular clouds, galactic rotation etc.) Only 
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Fig. 10. Left: Average surface density of star formation rate >I/ (in M© yr~^ pc~^) as 
a function of average gas (HI+H2) surface density Eg (in M© pc~^, in spirals (circles) 
and starbursts (squares); the solid line corresponds to ^' cx E^ **. Right: Average surface 
density of star formation rate 'if (in M0 yr~^ pc~^) as a function of ^a/Tdyn, where the 
dynamical timescal Tdyn ~ R/V (for rotational velocity V at radius R); the solid line 
corresponds to ^' oc 'Eg/ Tdyn (from Kennicut 1998). 



in the case of spirals predictions can be made, based on various instability criteria 
(e.g. Elmegreen 2002 for a review). In the case of other galaxies, a parametrised 
formulation of the SFR is adopted. Schmidt (1959) suggested that the density of 
the SFR is proportional to some power of the density of gas mass tog: 

•9 = ly (2.26) 

a formulation which has the merit of reminding us that stars are formed from gas, 
after all. However, it is not clear whether volume density p or surface density 
E should be used in Eq. 2.26. When comparing data with GCE models for the 
solar neighborhood, Schmidt (1959) uses surface densities (S in Mo/pc^). But, 
when finding "direct evidence for the value of TV" in his papei@ he uses volume 
densities {p in Mo/pc^) and finds N=2. Obviously, since S=/ p{z)dz, one has: 

j:^ ^J^piz^dz. 

It is not clear then which density should be used in the Schmidt SFR law: 
Volume density is more "physical" (denser regions collapse more easily) but surface 
density is more easily measured in galaxies. Furthermore, at first sight, it seems 



^Schmidt(1959) describes the distributions of gas and young stars perpendicularly to the 
galactic plane {z direction) in terms of volume densities poas oc exp{—z/hGas) and pstars oc 
exp(— z / h star s) with corresponding scaleheights (obervationally derived) hQas=7S pc and 
'iStars=144 pc~2 hcasi from that, Schmidt deduces that pstars <x Pq^s' f^^t is N=2. 
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that the density of molecular gas should be used (since stars are formed from 
molecular gas), and not the total gas density. 

Surprisingly enough, Kennicutt (1998) finds that, in normal spirals, the surface 
density of SFR correlates with atomic rather than with molecular gas; this con- 
clusion is based on average surface densities, i.e. the total SFR and gas amounts 
of a galaxy are divided by the corresponding surface aerea of the disk. In fact, 
Kennicut (1998) finds that a fairly good correlation exists between SFR and to- 
tal (i.e. atomic + molecular) gas. This correlations extends over four orders of 
magnitude in average gas surface density ps and over six orders of magnitude in 
average SFR surface density from normal spirals to active galactic nuclei and 
starbursts (see Fig. 10, left) and can be described as: 

* oc E^" (2.27) 

i.e. N=1.4. However, Kennicutt (1998) notes that the same data can be fitted 
equally well by a different N value, this time involving the dynamical timescale 
Tdyn = R/V{R)i where V{R) is the orbital velocity of the galaxy at the optical 
radius R: 

* (X (2.28) 

More recent observations confirm Kennicutt's results, but are also compatible with 
other SFR laws suggested in the literature (see e.g. Boissier et al. 2003). Note 
that, for a given data set, the determination of coefficient u in Eq. 2.26 depends 
on the determination of exponent N , and vice versa. 

In some cases, it is useful to consider the efficiency of star formation e, i.e. 
the SFR per unit mass of gas. In the case of a Schmidt law with A''=l one has: 
e = V = const., whereas in the case of N=2 one has: e = vmQ. 

Finally, one may ask whether the SFR law in Fig. 10 extends down to lower 
surface densities or whether there is a threshold of star formation in galactic disks. 
The data of Kennicutt (1998 and Fig. 10) apparently suggest a lower threshold 
of a few Mo/pc^. However, more recent data, obtained with the GALEX UV 
satellite, reveal star formation activity even in regions of gas density as low as 0.1 
Mq/pc^ (Fig- 11)- Those recent data imply that the threshold of star formation 
is much lower than thought before (and, perhaps, non existant). 

2.6 Gaseous flows and stellar feedback 

A galaxy is clearly not an isolated system, and it is expected to exchange matter 
(and energy) with its environment. This is true even for galaxies which are found 
away from galaxy groups. First of all, most of baryonic matter in the Universe 
today (and in past epochs) is in the form of gas residing in the intergalactic medium 
(e.g. Fukugita and Peebles 2004) and part of it is slowly accreted by galaxies. Also, 
small galaxies are often found in the tidal field of larger ones, and their tidal debris 
(gas and/or stars) may be captured by the latter. In both cases, gaseous matter 
is accreted by galaxies, and in the framework of the simple GCE model this is 
generically called infall. 
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Fig. 11. Surface density of star formation rate ^ as a function of gas (HI+H2) surface 
density Eg in galaxies observed with GALEX. The vertical band corresponds to the SFR 
treshold suggested by Kennicutt (1998). Clearly, UV observations show SFR activity 
below that threshold (from Boissier et al. 2006). 



On the other hand, gas may leave the galaxy, if it gets sufficient (1) kinetic 
energy or (2) thermal energy and (3) its velocity becomes larger than the escape 
velocity. Condition (1) may be met in the case of tidal stripping of gas in the field 
of a neighbor galaxy or in the case of ram pressure from the intergalactic medium. 
Condition (2) is provided by heating of the intestellar gas from the energy of 
supernova explosions, especialUy if collective effects (i.e. a large number of SN in 
a small volume, leading to a superbubble) become important. Finally, condition 
(3) is more easily met in the case of small galaxies, with swallow potential wells. 
Note that, since galaxies (i.e. baryons) are embedded in extended dark matter (non 
baryonic) haloes, a distinction should be made between gas leaving the galaxy but 
still remaining trapped in the dark halo and gas leaving even the dark halo. In 
the former case, gas may return back to the galaxy after "floating" for some time 
in the dark halo and suffering sufficient cooling. In the framework of the simple 
GCE model, all those cases are described generically as outflows. 

The rate of infall or outflow is difficult to calculate from first principles. In the 
case of infall, this is possible, in principle, for a hydrodynamical model evolving 
in an appropriate cosmological framework. In the case of outflows, the interac- 
tion between stars (SN) and the ISM, known as feedback, also requires detailed 
hydrodynamic modelling. No satisfactory models exist up to now for such com- 
plex processes. Incidentally, the treatment of feedback also affects the SFR of the 
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system (by making gas unavailable for star formation, either by heating it or by 
pushing it out of the system alltogether). 

In simple GCE models, infall and outflow are treated as free parameters, ad- 
justed as to reproduce observed features of the galaxian systems under study. 
Such features are the metallicity distributions (MD) of long -lived stars, or the 
mass-metallicity relationship of external galaxies. The former feature is nicely il- 
lustrated in the cases of the MD in the local disk and the halo (Sec. 3.1 and 3.2, 
respectively), which provide strong constraints on the history of those systems. 

2.7 Analytical solutions: the Instantaneous Recycling Approximation (IRA) 

The system of GCE equations (2.1 to 2.9) can be solved analytically if one adopts 

the Instantaneous Recycling Approximation (IRA), introduced by Schmidt (1963). 
Stars are divided in: "eternal" ones (those of low mass, with lifetimes far exceeding 
the age of the system) and "dead at birth" ones (massive stars, with lifetimes far 
shorter than the age of the system) for which it is assumed that tm=0. The 
dividing line between the two classes depends on the system's age, and for most 
practical purposes it is put at 1 Mq, corresponding to an age T ~12 Gyr. 

Assuming IRA allows one to replace ^'(t — tm) in the equations of GCE by 'l'(t) 
and thus to take out the SFR \l/(f) from the mass integrals. The equations can 
then be solved analytically (see e.g. Tinsley 1980). Solutions involve the Return 
fraction R defined in Eq. 2.25 and the yield pi of a given nuclide, defined as: 



The yield pi is the newly created amount of nuclide i by a stellar generation, per 
unit mass of stars blocked into "eternal" objects: indeed, yi{M) are the net yields 
of nuclide i (see Sec. 2.3.1) and 1 — i? is the mass blocked in low mass stars and 
compact objects (for the normalized IMF of Eq. 2.24). Obviously, the stellar yields 
Yi{M) and net yields yi{M) are properties of individual stars, while the yield pi 
is an integrated property of the IMF, as is the return fraction R. Note that only 
Yi{M) are interesting for numerical models, whereas yi{M) and Pi are useful only 
in analytical models. 

With those definitions, the gas mass fraction Xi of nuclide i in the case of the 
Closed Box model is: 



where a — nriQ/m is the gas fraction and A;_o the initial abmidance (X/.o=0 for 
metals) . This is the main result of IRA, relating the chemical enrichment of the 
gas to the amount of gas left. It is independent on time or on the form of the 
SFR, and for those reasons it is a powerful tool to study gas flows in the system. 
It can also be used to derive the metallicity distribution of stars and, thereof, the 
past history of the system (see Sec. 3.1.1 and 3.2.2). 




(2.29) 



Pi ln{ ) = Pi ln{—) 



(2.30) 
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Time (Gyr) Gas fraction (nig/m) 

Fig. 12. Results of calculations for a Closed Box model with Instantaneous Recycling 
Approximation (IRA, dotted curves) vs Non-IRA {solid curves). Left: Metallicity (upper 
curves) and gas fraction (lower curves) as a function of time. Right: Metallicity as a 
function of gas fraction. It is assumed that 5* = 1.2 mc Gyr"'^. 



If a simple Schmidt law of the form "if = vmc is adopted for the SFR (Eq. 
2.26) one can obtain the following solutions for the evolution of gas 

Tog = me-"'^--")* (2.31) 

and for the abundances 

X, ~ X,,o = P^v {l~R)t (2.32) 

which also satisfy the more general solution of Eq. 2.30. Thus, metallicity is 
roughly proportional to time, a result which is approximately valid even when 
IRA is relaxed. This property allows one to use stellar metallicity (especially iron, 
which has many strong and easily identifiable spectral lines) as a proxy for time, 
since stellar ages are notoriously difficult to evaluate. 

IRA turns out to be a surprisingly good aproximation for nuclides produced 
in massive stars, like e.g. oxygen, provided gas fractions stay above ~10%; this is 
illustrated in Fig. 12 for the case of a Closed Box model, with analytical solutions 
given by 2.31 and 2.32. Analytical solutions assuming IRA can also be obtained 
in the cases of gaseous flows into or out of the system. However, some specific 
assumptions have to be made about the form of those flows, and this limits a lot 
their interest. For instance, in the case of outflow at a rate proportional to the 
SFR, o = fc 5*, metalicity evolves as: 

X, - X,fl = ln{-) (2.33) 

1 + A; a 
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i.e. at the same gas fraction mctallicity is smaller than in the Closed box or, 
equivalently, a larger fraction of the system has to turn into stars in order to reach 
the same metallicity. Eq. 2.33 is formally the same as 2.30, with an effective yield 

Pi,eff = (2.34) 

in the place of the true yield pi, and Pi.eff < Pi- It is shown (Edmunds 1990) that 
gas flows always produce reduced effective yields, i.e. that metallicity increases 
most efliciently in the Closed Box model. For other analytical solutions, especially 
for models with specific forms of infall, see e.g. Pagel (1997) or Matteucci (2001). 

In the case of a secondary nuclide (see Sec. 2.3.2), its yield ps is proportional 
to the abundance of a primary one, i.e. ps = a Xp, and its evolution in the 
framework of Closed Box plus IRA is described by: 

1 a 

Xs = aXpln{-) = — X% (2.35) 
u pp 

that is, its abundance increases faster than the abundance of the primary (pp 
being the yield of the primary). 

The abundance ratio of two primaries or two secondaries remains constant 
during galactic evolution and equal to the ratio of the corresponding yields, in- 
dependently of the IMF and the SFR. Observations of that ratio allow then one 
to test theories of stellar nucleosynthesis, independently of the underlying GCE 
model. This is true only if both nuclides are produced in the same site (e.g. silicon 
and calcium, products of massive stars). On the other hand, an element produced 
mainly in a long-lived site may appear late in the Galaxy, even if it is a primary 
one. This is, among others, the case of Fe: its observed abundance increases faster 
than the one of oxygen after the first Gyr, and this increase is attributed not to a 
secondary nature but to a delayed production in SNIa (see Sec. 3.1.2). 

3 The Chemical evolution of the Solar Cylinder 

Obviously, the larger the number of observables, the more the framework pro- 
vided by GCE can be useful in constraining the history of the system. The Milky 
Way and, in particular, the Solar cylinder, is the best observed system today (for 
obvious reasons). 

The Solar cylinder may be defined as a cylindrical region of radius 0.5 kpc, 
perpendicular to the galactic plane and centered at Sun's position (at 8 kpc from 
the Galactic center, see Fig. 13). The interstellar gas is located near the plane and 
it is (usually assumed to be) chemically homogeneous. On the other hand, three 
stellar populations coexist, to various extents, in that region; they are distinguished 
by their kinematic and chemical properties. 

• The galactic halo: old (age ^12-13 Gyr), "pressure" supported, metal poor 
([Fe/H]<-1) with little or no rotation, large velocity dispersion perpendicu- 
larly to the disk and small contribution (<5%) to the overall surface density. 
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• The thin disk: young (age ~5-6 Gyr on average), rotationally supported, 
metal rich ([Fe/H]~-0.1 on average), with small scaleheight {h ^300 pc for 
the old stars, but only h ~100 pc for the young stars and the gaseous layer), 
dominating the total mass (75%). 

• The thick disk: older than the thin disk (^10 Gyr), more extended (scale- 
height h r^lA kpc), rotating, with moderate velocity dispersion, moderately 
metal poor ([Fe/H]'^-0.7) and contributing ~20% to the total surface density. 

It is tempting to assume, in the framework of the old, monolithic collapse scenario 
for the Milky Way's evolution (Eggen et al. 1962), a temporal continuity in the 
formation of those three components, i.e. halo — > thick disk — > thin disH^°l. Such 
a continuity is more difficult to establish in the modern framework of hierarchical 
galaxy formation: the three components may be totally uncorrelated, e.g. the 
thick disk may have been formed from tidal debris of satellite galaxies, while the 
thin disk from a slow accretion processed. In the following we shall study the 
cases of the local thin disk (often called the Solar Neighborhood) and of the Milky 
Way's halo independently. 



3.1 The local thin disk 

In the case of the Solar neighborhood, the number of available observational data 
is larger than for any other galactic system, allowing one to constrain strongly 
(albeit not in a unique way) the history of the system. Those data are (see e.g. 
Boissier and Prantzos 1999): 

• The current surface densities of gas (Eg ^^12 M© pc~^), live stars ((E* ^^30 
Mq pc~^), stellar residues (Sc ~8 Mq pc~^) and total amount of baryonic 
matter (E^ ~50 M© pc~^), as well as the current star formation rate (^o=2- 
5 Mq pc~^ Gyr~^); the corresponding gas fraction is a ^^0.24. 

• The elemental and isotopic abundances at solar birth {Xi^Q, i.e. it is as- 
sumed that the Sun's composition is typical of the ISM 4.5 Gyr ago) and 
today {Xi^Now)- Those two compositions are quite similar, suggesting little 
chemical evolution in the past 4.5 Gyr. 

• The stellar age-metallicity relationship (AMR), traced by the Fe abundance 
of long-lived stars, [Fe/H]=/(t) or log{Z) = f{t). 

• The metallicity distribution (MD) of long-lived G-type stars ^[j^J/^] or j^i^fj^^) 
showing that few of them were formed at [Fe/II]<-0.7 (0.2 Xpe.Q). 



^"In fact, it is difficult to account for the chemical and kinematical properties of all three 
galactic components in the monoUthic collapse scenario. 

^^The monolithic collapse scenario may describe better the properties of the inner halo and of 
the bulge. 
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Fig. 13. The various components of the Milky Way (bulge, halo, disk) and their main 
features; a distinction should be mde between thin and thick disk (not appearing in the 
figure). The Solar cylinder is at 8 kpc from the center (adapted from Pagel 1997). 



• The oxygen vs. iron (0-Fe) relationship, interpreted in terms of a delayed 
(after ~1 Gyr) enrichment of the ISM with products of SNIa, contributing 
to more than half of the solar Fe. 

Among those constraints, the age-metallicity relationship and the G-dwarf 
metallicity distribution are the most important ones. In principle, by combin- 
ing them, one may derive straightforwardly the SFR history [dn/dt) of the solar 
neighborhood through: 

dn dn d{logZ) 

dt d(logZ) dt 

In practice, however, this is impossible, because of the sensitivity of the result to 
the slope of the adopted AMR: a small variation in the form of the AMR produces 
a dramatic effect on the resulting SFR history (see Prantzos 1997). For that 
reason, the local SFR history is reconstructed only indirectly: models of the GCE 
of the Solar neighborhood are developed, which must satisfy all the aformentioned 
observational constraints. Among those constraints, the metallicity distribution is 
probably the most significant. 
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3.1.1 The Local Metallicity Distribution of long-lived stars 

In the framework of the Closed box model of GCE with IRA, the MD of long-lived 
stars can be derived as follows. In order to reach a given metallicity Z , a certain 
amount of stars has to be created, given by (Eq. 2.4): ms — m — ma, or (by 
normalising to the total mass m): 

n = 1 - <T (3.2) 

where a — mG/m is the gas fraction and n = ms/m the star fraction. For a 
system with a final metallicity Zi and star fraction ni — 1 — (Ti, the cumulative 
metallicity distribution (CMD), i.e. the number of stars with metallicity lower 
than Z as a function of Z, is given by 



n{< Z) 1 



(3.3) 



By using the fundamental result of IRA, namely a = exp{—Z/p) (Eq. 2.30) 
and taking the derivatives of Eq. 3.3, one obtains the differential metallicity 
distribution (DMD): 



din/ni) _ ln{lO) Z 

d(logZ) 1 — exp{—Zi/p) p 



- e-^/P (3.4) 



i.e. the number of stars per logarithmic metallicity interval as a function of metal- 
licity Z. This relation has a maximum for Z = p, i.e. when the metallicity is 
equal to the yield (both metallicity and yield can be expressed in units of the solar 
abundance of the corresponding element). 

Eq. 3.4 for the Closed box model with initial metallicity Zq=0 appears in Fig. 
14, where it is compared to data for the local disk. It predicts many more stars at 
low metallicities than observed, a problem known as the "G-dwarf problem'f^. 

Two of the main solutions proposed for the G-dwarf problem appear also on Fig. 
14. According to the first one, the disk started with an initial metallicity Zq ^0.1 
ZQ^pre-enrichment) . In that case, all metallicities in Eq. 3.4 are replaced by Z—Zq 
and the resulting curve fits relatively well the data. The main drawback of that 
hypothesis is that it is hard to justify the origin of such a large pre-enrichment: it 
is true that the Galactic halo, which preceded disk formation, reached a maximum 
metallicity of ~0.1 Zq (for Fe); but its average (stellar) metallicity is ~0.03 Z0(for 
Fe) amd its total mass ('^2 10^ M©) is 20 times smaller than the one of the disk 
(~4.5 10^° Mg). There is simply not enough mass and metals produced in the 
halo to justify pre-enrichment of the disk to such a high level. Moreover, the halo 
has a low specific angular momentum (contrary to the disk) and material escaping 
it should be accreted rather by the bulge, not the disk. 



G-type stars are bright enough for a reasonably complete sample to be constructed and long- 
lived enough to survive since the earliest days of the disk; the same problem is encountered if F- 
or K- type stars are used. 
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Fig. 14. Metallicity distribution of long-lived stars for three models: Closed Box, ex- 
ponentially decreasing Infall (with a timescale of 7 Gyr) and Pre-enrichment (with 
Xo ~ O.O8X0 for Fe), respectively. They are compared to data for the Solar neigh- 
borhood. 



The second hypothesis is that the disk did not evolve as a Closed box, but was 
gradually built from infall of metal free (or metal poor) material. In the Closed 
box, all the gas of the system is available from the very beginning; a large stellar 
activity is then required to enrich all that gas to, say, 0.1 Zq, and correspondingly 
many small and long-lived stars are formed at low Z. In the case of infall, only a 
small amount of gas exists early on; it takes then a small number of SN to enrich 
it to 0.1 Z0, and correspondingly few low mass stars are formed at low Z. 

Infall appears then as an elegant solution to the local G-dwarf problem, es- 
pecially in view of the fact that gas accretion to galaxies is expected to be a 
common phenomenon in the Universe. The rate of the infall is not precisely de- 
termined by the data of the local disk. An exponentially decreasing infall rate 
f{t) = A exp{—t/T) with a long characteristic time-scale of r r^7 Gyr provides 
a reasonably good fit to the data and, in view of its simplicity, is often used in 
models of the solar neighborhood (see next section). However, other forms may 
do as well and even better; this is the case, for instance, of a gaussian (as a func- 
tion of time) infall rate, with a maximum prior to solar system formation (see e.g. 
Prantzos and Silk 1998). 

The infall rate f(t) must obviously be normalized such as: 




(3.5) 



where T,t is the total surface density in the Solar neighborhood (see Sec. 3.1) and 
T ~12 Gyr is the age of the Galactic disk. 
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Fig. 15. History of the Solar neighborhood, according to a GCE model with infall, 
constrained by various observables. Left: Results are plotted as a function of time (or 
age, for the bottom panel). Data for the present-day local disk are displayed with vertical 
bars. Data for the age-metallicity relation are from various sources. Right: Results are 
plotted as a function of metallicity [Fe/H]. See text for comments on the various curves. 



3.1.2 A brief History of the Solar neighborhood 

The observed properties of the local disk, presented in Sec. 3.1, along with the 
constraint of G-dwarf MD (analysed in Sec. 3.1.1), "dictate" the parameters of 
simple GCE models that may be built for that system. The results of such a model 
are displayed in Fig. 15, as a function of time (left panels) and of metallicity [Fe/H] 
(right panels). The various paremeters of the model are adjusted as follows: 

• The total amount of infalling matter is normalized to the local surface density 
through Eq. 3.5. 

• The time-scale of the adopted exponentially decreasing infall rate is suffi- 
ciently long (t= 7 Gyr) as to reproduce the metallicity distribution. 



The coefficient i/=0.3 Gyr^^ of the SFR * = mc (Eq. 2.26) is adjusted 
as to leave the system at T—12 Gyr with a gas fraction a ^0.2, as observed. 
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Fig. 16. Abundances of the chemical evolution model of the solar neighborhood 4.5 Gyr 
ago, compared to solar abundances, from C to Zn. Elemental abundances (top) are co- 
produced within a factor of 2. Isotopic abundances (bottom) are also well reproduced, 
except ^^N (for which novae -not included in this calculation - appear to be the source) 
and some Ca and Sc isotopes. Note that mono-isotopic F is produced by neutrino- 
nucleosynthesis in massive stars, according to the calculations of Woosley and Weaver 
(1995) adopted here. More than half of solar Fe results from SNIa (with yields adopted 
form Iwamoto et al. 1999), while part of C and N comes from intermediate mass stars 
(with yields adopted from van den Hoek and Gronewegen 1997). 



With those parameter adjustments, it remains to be seen whether the other 
observables of the system (Sec. 3.1) are reproduced. In Fig. 15 it is seen that 
the current model SFR is well within observational uncertainties. Note that the 
derived SFR history is rather flat, around an average value of ^3.8 Mq/jt (alter- 
natively, it could be fitted with a broad gaussian). 

The resulting age-metallicity relationship (left bottom panel in Fig. 15) fits the 
data approximately, but its should be stressed that uncertainties in stellar ages are 
fairly large; in fact, dispersion in the age-metallicity relation is even larger than 
displayed in Fig. 15, to the point that this observable may be of very little use as 
a constraint for the local GCE (see Sec. 3.1.4 and Fig. 18). 
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The rise of Fe/0 (right panel of Fig. 15) is due to delayed contribution of 
Fe by SNIa. The adopted SNIa rate is from Greggio and Renzini (1983) and it 
is assumed that 0.04 of the binaries produce SNI; the resulting SNIa rate as a 
function of time appears on the left panel. 

Finally, the local GCE model, combined with the adopted stellar yields, should 
also reproduce the pre-solar composition, well established after meteoritic and 
photospheric measurements (e.g. Loddcrs 2003). The results of such a comparison 
(model results at time r-4.5 Gyr divided by pre-solar composition) are displayed 
in Fig. 16. It can be seen that all elements and almost all isotopes are nicely 
co-produced (within a factor of 2 from their pre-solar value), with key elements 
such as O and Fe being very well reproduced. Taking into account the large 
abundance variations between O and Sc (a factor of 10^) this agreement should 
be considered as a triumph for stellar nucleosynthesis models. Some exceptions to 
this achievement are noted in the legend of Fig. 15, but the overall result is quite 
encouraging: the adopted IMF, SFR and stellar yields reproduce very well the 
solar system composition (assumed to be typical of the local ISM 4.5 Gyr ago). 

3.1.3 The local evolution of Deuterium 

Modelling the GCE of deuterium is a most straightforward enterprise, since this 
fragile isotope is 100% destroyed in stars of all masses (already on the pre-main 
sequence) and has no known source of substantial production other than Big Bang 
nucleosynthesis (BBN). In the framework of the simple Closed box model with 
IRA, its abundance is given by: 

D = Dp a-T^ (3.6) 

where Dp is the primordial abundance, a the gas fraction and R the return frac- 
tion. If the boundary conditions of its evolution (namely the primordial abundance 
Dp and the present day one) were precisely known, the degree of astration, i.e the 
degree of processing of galactic gas inside stars, which depends on the adopted IMF 
and SFR, would be severely constrained. This, in turn, would be another strong 
constraint (not appearing in the list of Sec. 3.1) for models of local GCE. For 
typical values of tr ^ 0.2 and R ^0.3 one obtains: D/Dp ~0.5 in the Closed box 
model, i.e. depletion by a factor of ~2. However, the true depletion factor should 
be smaller, because infall (of presumably primordial composition) is required for 
the GCE of the local disk, as argued in the previous section. 

The primordial abundance of D is now well determined, since observations of D 
in high redshift gas clouds agree with abundances derived from observations of the 
Cosmic Microwave Background combined to calculations of standard BBN (see e.g. 
Prantzos 2007 and rc;fcrcncc therein). However, the present day abundance of D in 
the local ISM is now under debate. Indeed, UV measurements of the FUSE satellite 
along various lines of sight suggest substantial differences (a factor of two to three) 
in D abundance between the Local Bubble and beyond it. Until the origin of that 
discrepancy is found (see Hebrard et al. 2005 and Linsky et al. 2006), the local 
GCE of D in the past few Gyr will remain poorly understood: naively, one may 
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Fig. 17. Evolution of deuterium in the solar neighborhood, as a function of time. The 
adopted model satisfies all major local observational constraints (see Fig. 15). Data 
(with error bars or inside box) correspond to: primordial Dp (resulting from standard 
Big Bang nucleosynthesis + baryonic density inferred from analysis of the cosmic mi- 
crowave background), observations of high redshift gas clouds (DLAs), pre-solar (inferred 
from observations of stellar wind composition) and in the local ISM (through different 
interpretations of UV data from the FUSE satellite). 

expect that a high value would imply strong late infall of primordial composition, 
while a low value would imply strong late astration. In any case, corresponding 
models should also satisfy all other local observables, like the overall metallicity 
evolution and the G-dwarf metallicity distribution, which is not an easy task. 

3.1.4 Uncertainties in the evolution of the Solar neighborhood 

Our understanding of the evolution of the solar neighborhood depends on how well 
are established the relevant observational constraints (listed in Sec. 3.1). It has 
recently been suggested that some of those constraints may be less well understood 
than thought before. 

Stellar ages are much harder to evaluate than stellar metallicities, and the 
form of the local AMR has varied considerably over the years. The seminal work 
of Edvardsson et al. (1993) on 189 F-dwarfs established a clear trend of decreasing 
metallicity with age, albeit with substantial scatter (Fig. 18 right). Such a trend is 
compatible with (and predicted by) all simple models of local GCE, either closed 
or open (i.e. with infall) models. It should be noted, however, that the adopted 
selection criteria in that paper introduced a bias against old metal-rich and young 
metal-poor stars. 
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Fig. 18. Age-metallicity relationship in the solar neighborhood. Left: data from Nord- 
strom et al. (2004) for FGK stars within 40 pc. Thick dots indicate average metallicities 
in the corresponding age bins; no age-metallicity relation appears in those data. Right: 
data for 189 F dwarfs from the study of Edvardsson et al. (1993), suggesting an age- 
metallicity relation; by construction, however, old and metal-rich stars are excluded from 
that sample. Figure from Nordstrom et al. (2004). 



The large survey of Nordstrom et al. (2004), concerning ^^14000 F and G stars 
with 3-D kinematic information (but less precise spectroscopy than the Edvardsson 
et al. study), provides a radically different picture: the volume limited subsample 
of 462 stars with " well-defined" ages withing 40 pc displays a flat AMR (an average 
metallicity of [Fe/H]~-0.2 at all ages) with a very large scatter (Fig. 18 left). 
Acounting for the fact that the oldest stars have the largest age uncertainties 
does not modify the flatness of the AMR. If confirmed, a flat AMR would require 
different assumptions than those adopted in current models (e.g. substantial late 
infall to dilute metals, with a strong impact on D evolution or on the metallicity 
distribution). 

The form of the second "pillar" of local GCE, namely the local MD, has also 
been revisited recently by Haywood (2006). He argues that, since the more metal- 
poor stellar populations have larger dispersion velocities vertically to the disk, 
larger scaleheight corrections should be applied to their numbers in order to get 
their true surface density. By adopting such steeply decreasing correction factors 
with metallicity (Fig. 19 left) he finds then that a closed box model nicely account 
for the corrected MD (Fig. 19 right). He notes that the (volume corrected) local 
population at -l<[Fe/H]<-0.5 comprises 15-20 % of the total; it should then be 
"naturally" considered as the thick disk population, which weights ~7 Mq pc~^, 
compared to ~35 Mq pc~^ for the total stellar population in the solar cyhnder. 

The relationship of the thick and thin disks is not yeat clear; if thick disk stars 
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Fig. 19. Left: Correcting factors as a function of metallicity, to be applied to the local 
(thin disk) metallicity distribution, according to Haywood (2006); they are based on 
estimates of velocity dispersion as a function of metallicity. Right: If (the inverse of) 
those correction factors are applied to the results of a closed box model, the resulting 
metallicity distribution (grey shaded aerea) ressembles closely the observed one (from 
Haywood 2006). 



were mostly accreted from merging/disrupted satellites, that population cannot 
be considered as belonging to the early phase of the disk and the arguments of 
Haywood (2006) do not hold. Also, in view of the large scaleheight of the thick 
disk (1400 pc) and of the elliptical orbits of its stars (vs. circular orbits for those 
of the thin disk) one may wonder how large the local "chemical box" could be 
and still be considered as a single system with well defined evolution. It is clear, 
however, that corrections to the observed local MD should be carefully considered 
before comparing to GCE models and this would certainly impact on the required 
infall timescales (which could be smaller than the "canonical" value of ~7 Gyr 
currently adopted). 

3.2 The Milky Way halo 

The list of observational constraints for the Galactic halo is much shorter than 
for the solar neighborhood. In particular, we have no information on the gas 
fraction at the end of its evolution, no isotopic abundances available and there is 
no age-mctallicity relationship. Age differences of Galactic globular clusters are 
evaluated to ~l-2 Gyr (e.g. Salaris and Weiss 1993) and this is usually taken to 
be the duration of halo formation. Note that this timescale is much longer than 
the ~10^ yr derived in the monolithic collapse scenario of Eggen et al. (1962). 

The two most important pieces of information on the Galactic halo are: (i) the 
abundance patterns of various elements and (ii) the metallicity distribution. 
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3.2.1 Halo star composition: any signatures of Pop. Ill stars ? 

Detailed abundance patterns have now been established in the case of halo stars, 
for the majority of the chemical elements (see Hill, this volume). In the case of 
heavier than Fe elements, their interpretation is hampered at present by our poor 
understanding of the production sites of the s- and r- nuclides (see Arnould, this 
volume). In the case of intermediate mass elements, up to the Fe-peak, the sit- 
uation is more favourable than in the case of the local disk, since the system is 
thought to be enriched only by the products of massive stars and not of SNIa 
(which operate effectively on longer timescales than ~1 Gyr). Observed abun- 
dance patterns of those elements constrain then corresponding massive star yields 
(folded with a stellar IMF). In particular, the composition of the most metal-poor 
stars (say, below [Fe/H]~-3) constrains the yields of the first massive star gener- 
ations (see also Limongi, this volume), those formed from material of primordial 
composition (without metals) and generically called Pop. Ill stars. 

On theoretical grounds, it has been conjectured that Pop. Ill stars should be 
more massive on average than stars of subsequent, metal enriched, generations. 
The reason is that, in the absence of metals, gas cannot radiate efficiently away 
its heat (the most efficient cooling agent in that case being the H2 molecule); 
thus, it cannot reduce sufficiently its internal pressure and collapse gravitationally. 
Even in the absenc;e of metals, the complexity of the problem of gas collapse and 
fragmentation (involving various forms of instabilities, turbulence, magnetic fields, 
ambipolar diffusion etc.) is such that the typical mass of the first stars is very 
poorly known at present (sec Glover 2005 and references therein). It may be as 
"low" as the mass of a typical massive O star today (i.e. a few tens of Mq) 
or as large as several hundreds of Mq. It is even possible that, depending on 
ambient density, two ranges of typical masses co-exist for Pop. Ill stars, i.e. one 
around a few M©, and another around a few hundreds of (Nakamura and 
Umemura 2002). Is there any evidence from current observations for or against 
those theoKrtieal ick^as ? 

A summary of observational data for intermediate mass elements in halo stars 
appears in Fig. 20, where comparison is made with a GCE model (from Goswami 
and Prantzos 2000). The adopted stellar yields (from Woosley and Weaver 1995 
and Chieffi and Limongi 2004) are from non-rotating Core Collapse supernovae 
(CCSN) in the "normal" mass range of 12 to ~40 M©, exploding with "canonical" 
energies of E =10'''^erg and with progenitor metallicities from Z=0 to Z=Z0. The 
salient features of that comparison can be summarized as follows: 

• The a/Fe ratio (where a stands for alpha elements like O, Mg, Si, Ca) 
remains ^constant^2-3 times solar from [Fe/H]^-1 down to the lowest ob- 
served metallicities. This is attributed to SNII, while its solar value is at- 
tributed to the later contribution of Fe by SNIa. 

• In the adopted stellar models (no mass loss or rotation) N is produced as 
secondary, while observations show that it behaves as a primary (i.e. N/F is 
constant); this point will be discussed in the next section. 
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Fig. 20. Abundance ratios [X/Fe] as a function of metallicity [Fe/H] in stars of the 
Milky Way; small data points are from various sources, while the large data points at 
low metallicity are from the VLT survey of Cayrel et al. (2004). Results of standard 
chemical evolution models, performed with two sets of metallicity dependent massive 
star yields (Woosley and Weaver 1995, solid curves ; Chieffi and Limongi 2004, dashed 
curves) are also displayed. Yields for SNIa are adopted from Iwamoto et al. (1999) 
and for intermediate mass stars from van den Hoek and Gronewegen (1997). While the 
behaviour of alpha elements and Mn is correctly reproduced (at least qualitatively), there 
are large discrepancies in the cases of Cr, Co and Zn. 



• Abundances of elements affected by the odd-even effect (see Sec. 2.3.2), like 
Na and Al, require uncertain and potentially large NLTE corrections and it 
is difficult, at present, to use them in order to constrain stellar yields. 

• There is a group of elements below the Fe-peak, including K, Sc, Ti and V, 
the abundance pattern of which is poorly described by the adopted yields (as 
well as by those of the Tokyo group of Nomoto et al. 2006); this apparently 
"generic" problem of current stellar models is discussed by Limongi (this 
volume). 
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Fig. 21. Nucleosynthesis yields of a 25 Mq star, compared to observations of extremely 
metal-poor stars (vertical error bars, from Cayrel et al. 2004). Left: The star explodes 
with a canonical energy of 1.5 10^^ erg; Right: The star explodes as a hypernova, with 
an energy of 10 lO''^ erg (both figures from Nomoto et al. 2006). Abundance ratios in 
the Fe-peak are much better reproduced in the second case. 



• The decreasing trend of Mn/Fe and Cu/Fe with decreasing metallicity is in 
nice quantitative agreement with theoretical expectations. 

• The Fe-group nuclei Cr, Co and Zn behave very differently than expected 
on theoretical grounds. 

The reasons for the latter discrepancy probably lie in our poor understanding 
of the explosion mechanism and of the nature of the early CCSN. A large amount 
of work was devoted in recent years to understanding the peculiar abundance 
patterns among Fe-peak elements observed in extremely metal-poor (EMP) stars 
(for reviews see Cayrel 2006 or Limongi, this volume, and references therein). In 
those works it was explored whether some specific property of the Pop. Ill stars 
(characteristic mass, rotation, explosion energy etc.) was quite different from the 
corresponding one of their more metal rich counterparts and produced a different 
nucleosynthetic pattern in the ejecta. Two main classes of solutions (not mutually 
exclusive) appear to arise from those works: 

1) Normal massive stars (10-50 Mq), exploding as core collapse SN (CCSN) 
and leaving behind a neutron star or a black hole. It was shown (essentially by the 
Tokyo group, see Nomoto et al. 2006 for a recent review) that higher energies than 
the canonical one of i? = 10^^ ergs combined with asphericity of the explosion can 
help to improve the situation concerning the Zn/Fe and Co/Fe ratios observed in 
EMPs (Fig. 21 right) . Material ejected along the rotation axis (in the form of jets) 
has high entropy and is found to be enriched in products of a-rich freeze-out (Zn 
and Co), as well as Sc (which is generically undeproduced in spherical models); 
this kind of models seems at present the most promising, but this is not quite 
unexpected (since they have at least one more degree of freedom w.r.t. spherical 
models). Such properties (high energies, asphericity) are indeed observed in the 
local Universe for hypernovae, a class of energetic CCSN. However, hypernovae are 
not so frequent today as to affect the later chemical evolution of galaxies (normal 
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Fig. 22. Nucleosynthesis yields of Very Massive Stars (VMS, with mass >100 M©). 
Left: Stars in the 140 to 260 M© range explode as Pair Instability Supernovae (PISN) 
and do not produce enough Zn (from Hcgcr and Wooslcy (2002). Right: Stars in the 
300-1000 M© range explode after core collapse and produce significant amounts of Zn 
(from Ohkubo et al. 2006). Abundance ratios in the Fe-peak are much better reproduced 
in the second case. A Pop. Ill population composed exclusively of VMS (in the 140-1000 
M© range) and with an appropriate IMF is obviously compatible with observations of 
EMP stars. 



CCSN plus SNIa can easily account for e.g. the solar abundance pattern). If 
hypernovae indeed affected the early chemical evolution, they must have been 
much more abundant in Pop. Ill than today. No satisfactory explanation for 
such a large hypernova fraction in the early Universe exists up to now. It seems 
hard to conceive that the typical energy of the explosion was substantially higher 
in Pop. Ill stars (all other parameters, except metallicity, being kept the same). 
On the other hand, rotation effects may indeed be stronger in low metallicity 
environments, because of smaller angular momentum losses due to lower mass loss 
rates (e.g. Meynet, this volume). In any case, none of the models explored so 
far appears to account for the early trend of Cr/Fe, which remains a mystery at 
present. 

2) Very massive stars (VMS), above 100 Mq. Such stars are thought to: 
collapse to black holes, if M<140 Mq; explode during oxygen burning as pair- 
implosion SN (PISN), if 14O<M/M0<3OO; and again collapse without ejecting 
metals if M>300 Mq (sc!c Hciger and Woosley 2002 and references therein). In the 
case of PISN, it was shown that they do not produce enough Zn to account for 
the oberved high Zn/Fe ratio in EMP stars (Heger and Woosley 2002, Umeda and 
Nomoto 2006 and Fig. 22 left) and they were thus excluded as major contributors 
to the early chemical evolution of the MW. This shortc;oming cast doubt as to the 
existence of a Pop. Ill generation composed exclusively of VMS. However, in a 
recent work Ohkubo et al. (2006) explore the conditions under which rotating stars 
of M=500 and 1000 Mq do explode and eject metals (although a massive black hole 
is ultimately produced in both cases). They find explosions for a certain region of 
their parameter space and substantial production of Zn as a generic feature (see 
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Fig. 23. Evolution of oxygen, nitrogen and carbon in the Milky Way (from top to 
bottom). Observations are from various sources, while model results correspond to three 
sets of the Geneva stellar yields: stars with no rotation, with rotational velocity V=300 
km/s (as observed in the solar neighborhood) at all metallicities, and for V=800 km/s 
for low metallicity stars. The latter case provides the required primary N to fit the low 
metallicity data (see text). Timescales corresponding to metallicity and masses of stars 
dying on such timescales appear in the bottom panel. 



Fig. 22 right). Even if each class of PISN or 500-1000 Mq stars cannot reproduce 
the observed EMP abundance pattern alone, it is obvious that considering the 
full 140-1000 Mq range (folded with an appropriate IMF) will account for such 
a pattern. Thus, VMS cannot be excluded at present as candidates for Pop. Ill 
stars, at least not on "chemical" grounds. 

3.2.2 The quest for primary Nitrogen 

The behaviour of N as primary (i.e. [N/Fe]^0) was known for sometime, but it 
was recently confirmed from VLT mesurements (Spite et al. 2005) down to the 
realm of very low metallicities ([Fe/II]~-3, see Fig. 23 middle). For a long time, 
the only known source of primary N was Hot Bottom Burning (HBB) in AGB stars 
(see Sec. 2.3.4). Even the most massive AGBs (~8 M©) have lifetimes (^5 10^ 
yr) considerably longer than those of typical SNII progenitors (20 Mq stars living 



42 



Stellar Nucleosynthesis: 50 years after B2FH 



for '^lO'' yr) and it is improbable that they contributed to the earliest enrichment 
of the Galaxjf^. On the other hand, massive stars were though to produce N only 
as secondary (from the initial CNO) and not to be at the origin of the observed 
behviour. 

Rotationally induced mixing in massive stars changed the situation consider- 
ably: N is now produced by H-burning of C and O produced inside the stoiP^. i.e. 
as primary. The first set of models, rotating at 300 km/s at all metallicities, did 
not provide enough primary N at low metallicities to explain the data (see Fig. 
23). Assuming that low metallicity massive stars were rotating faster than their 
high-metallicity present-day counterparts (at 800 km/s) leads to a large produc- 
tion of primary N, even at low Z (see Meynet, this volume) and allows one to 
explain the data (Fig. 23). Thus, there appears to be a "natural" solution to the 
problem of early primary N, which may impact on other isotopes as well (e.g. ^^C, 
produced in a similar way). 

3.2.3 The early MW and hierarchical galaxy formation 

In the previous sections, the early chemical evolution of the MW was discussed 
independently of the cosmological framework in which it took place. According 
to the currently dominant paradigm of hierarchical structure formation, the early 
phases of a galaxy's evolution are the most complex ones, as they involve multiple 
mergers of smaller sub-units. In the case of the Milky Way, interesting "chemical 
signatures" of that period should still be left around us today, in the form of the 
metallicity distribution (MD) of long-lived stars . 

The MD of Galactic halo field stars (HMD) is rather well known in the metal- 
licity range -2.2<[Fe/H]<-0.8, while at lower metallicities its precise form has still 
to be established by ongoing surveys (Fig. 24 left). Its overall shape is well fitted 
by a simple model of GCE as dn/dlogZ cx Z/p e~^^P, where p is the yield of a 
stellar generation (Eq. 3.4); this function has a maximum for Z — p. The HMD 
peaks at a metallicity [Fe/H]=-1.6 (or [0/H]=-l.l, assuming [O/Fe]~0.5 for halo 
stars), pointing to a low yield p=l/13 of the corresponding value for the solar 
neighborhood. Such a reducedhalo yield (see Eq. 2.34) is "classicaly" (i.e. in the 
monolithic collapse scenario) interpreted as due to outflow during halo formation 
(Hartwick 1976), at the large rate of 8 times the SFR (Prantzos 2003). How can 
it be understood in the framework of hierarchical merging ? 

It should be noted that the typical halo metallicity ([Fe/H]=-1.6) is substan- 
tially lower, by more than 0.5 dcx, than the corresponding metallicities of nearby 
galaxies of similar mass [Mnaio ~2 10^ M©), as can be seen in Fig. 24 (right). 
That figure also displays the well known galaxian relationship between stellar mass 



^''Note, however, that the timescales of the early Galactic evolution are not constrained (there 
is no age-metalicity relation) and the contribution of AGBs even as early as [Fo/H]~-3 cannot 
be absolutely excluded. 

is produced after mixing of protons in He-rich zones, where ^^C is produced from the 3-a 
reaction. 
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Fig. 24. Left: Metallicity distribution of field halo stars from Ryan and Norris (1991, 
dots), and the ongoing research of Beers and collaborators (T. Beers, private communi- 
cation, asterisks). The curve is a simple model with outflow rate equal to 8 times the star 
formation rate. Right: Stellar metallicity vs stellar mass for nearby galaxies; data and 
model are from Dekel and Woo (2003). The MW halo, with average metallicity[Fe/H]=- 
1.6 (see left panel) or [0/H]=-l.l and estimated mass 2 10^ M© falls below that rela- 
tionship. 



and stellar metallicity; most probably, that relationship results from mass loss, 
which is more important in lower mass galaxies, since the hot supernova ejecta 
escape more easily their swallow potential well (e.g. Dekel and Woo 2003). 

Assuming that the MW halo has been assembled from sub-units similar to the 
low mass galaxies of Fig. 24 (right), one may interpret the HMD as the sum of 
the MD of such low mass galaxies; it is assumed that each one of those galaxies 
(of mass to) evolved with an appropriate outflow rate and corresponding effective 
yield p{m) = {1 — R)/{1 + k — R), where k{m) is the outflow rate in units of the 
SFR and R is the return mass fraction, depending on the adopted stellar IMF. 
In that case, one has: HMD(Z) = l/mnaio J Z/p{m) e-^/P'^'^^ (l){m)mdm, where 
fpim) is the mass function of the sub-units and p{m) is the effective yield of each 
sub-unit (obtained directly from Fig. 24 (right) as p{ni) = Z{m), i.e. smaller 
galaxies suffered heavier mass loss). 

The results of such a simple toy-model for the HMD appear in Fig. 25 (top 
and middle left, from Prantzos 2006). The HMD is extremely well reproduced, 
down to the lowest metallicities, assuming (j){m) cx m~^; such a halo mass function 
results from recent high resolution numerical simulations for Milky Way sized dark 
haloes (Diemand et al. 2006, Salvador! et al. 2006). Low mass satellites (down 
to 10^ Mq) contribute most of the low metallicity stars of the MW halo, whereas 
the high metallicity stars originate in a couple of massive satellites with M ~ 10^ 
Mq (Fig. 25 left, bottom). 
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Fig. 25. Left, top and middle panels: Metallicity distribution (in lin and log scales, 
respectively) of the MW halo, assumed to be composed of a population of smaller units 
(sub-haloes). The individual MDs of a few sub-haloes, from 10^ Mq to 2 10* M0, are 
indicated in both panels (but clearly seen only in the middle), as well as the sum over 
all haloes (upper curves in both panels, compared to observations). Small sub-haloes 
contribute the largest fraction of the lowest metalicity stars {bottom left). Right panels: 
Properties of the sub-haloes as a function of their mass. 



Some properties of the sub-haloes as a function of their mass appear also in 
Fig. 25 (right panels). The outflow rate, in units of the corresponding SFR, is 
k{m) = {1— R){pTrue I pirn) — 1), where pTrue is the yield in the solar neighborhood 
{pTrue = 0.7^0, from the local G-dwarf MD). If the MW halo were formed in a a 
potential well as deep as those of comparable mass galaxies, then the large outflow 
rate required to justify the HMD {k=8) is puzzling; on the contrary, the HMD 
is readily understood if the MW halo is formed from a large number of smaller 
satellites, each one of them having suffered heavy mass loss according to the simple 
outflow model. 

In summary, the halo MD can be understood in terms of hierarchical galaxy 
formation, provided that each merging sub-halo behaves as described by the simple 
GGE model with outflow^ the outflow rate depending on the sub-halo mass as 
suggested by observations of local dwarf galaxies. 
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4 Conclusion 

In this tutorial, we present the formalism and the ingredients of standard GCE 
models. We also illustrate the concepts of GCE by applying simple models to the 
study of the local disk and of the Milky Way's halo, two systems with a large 
body of available observational data. In both cases, useful information on those 
galaxian systems may be extracted from the data, interpreted in the framework of 
such simple models. It is stressed, however, that the history of the systems cannot 
be reconstructed in a unique way, because (i) dispersion in the data does not allow 
it, and/or (ii) some of the data (e.g. the metallicity distributions) are independent 
of the form of the star formation rate. 

It is clear that the future of the discipline lies in forthcoming chemodynamical 
models of GCE, incorporating several physical ingredients and evolving in a cosmo- 
logical framework. However, untill such models are properly validated/calibrated, 
simple GCE models such as those presented here will continue to be useful, allow- 
ing for a first, rapid interpretation of a large amount of observations. 
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